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PREFACE 


This  report  contains  an  invited  paper  given  at  the  Conference  on 
Natural  Physical  Sources  of  Underwater  Sound  at  the  University  of 
Cambridge,  July  1990.  It  was  written  under  the  partial  sponsorship  of 
Dr.  M.  Orr,  Office  of  Naval  Research  (Code  11250A),  and  Dr.  K.  Lima, 
Naval  Underwater  Systems  Center  (Code  10). 

The  results  presented  in  this  paper  are  the  cummulation  of  work 
begun  in  1980  and  supplement  the  results  presented  at  the  Nato 
Advanced  Workshop  on  the  Natural  Mechanisms  of  Surface  Generated 
Noise  in  the  Ocean  held  in  Lerici,  Italy,  in  1987,  and  published  in  Sea 
Surface  Sound.  The  paper,  "Low  Frequency  Ocean  Ambient  Noise: 
Measurement  and  Theory,"  is  included  here  as  an  appendix. 
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LOW  FREQUENCY  NOISE  FROM  BREAKING  WAVES 


William  M.  Carey  and  James  W.  Fitzgerald* 
Naval  Underwater  Systems  Center 
New  London,  Connecticut  06320 
USA 


ABSTRACT.  Recent  experiments  confirm  the  production  of  sound  by 
breaking  waves  at  lower  frequencies  (30  to  500  Hz).  Individual  breakers 
produce  impact  noise  as  well  as  a  random  collection  of  individual 
spectral  events.  Measured  ocean  ambient  noise  spectrum  levels  increase 
at  less  than  1  dB  per  octave  toward  a  broad  maximum,  which  has  a  weak 
wind  speed  dependence  between  300  to  500  Hz.  Noise  intensities  (<  500 
Hz)  are  a  function  of  wind  speed  (U)  to  the  2n  power  with  1.3  <  n  <  2.5 
and  a  value  of  n=1.5  at  200  Hz.  The  production  of  noise  in  this  region 
has  a  dipole  characteristic.  Breaking  waves  produce  an  impact,  bubble 
plume,  and  bubble  cloud.  The  dynamic  evolution  of  these  plumes  and 
clouds  provides  a  mechanism  for  sound  production.  Since  the  initial 
plume  and  cloud  have  appreciable  void  fractions,  compressible  resonant 
behavior  of  these  structures  as  a  whole  or  as  multiply  connected 
regions  can  be  represented  as  compact  acoustic  monopoles  and  dipoles. 
The  pressure  release  surface  would  result  in  an  effective  dipole 
characteristic.  Sufficient  energy  exists  in  the  initial  breaking 
vorticity  and  turbulence  to  explain  measured  source  levels.  Since  a 
good  radiator  of  sound  is  also  a  scatterer  of  sound,  these  plumes  and 
clouds  will  also  scatter  sound. 


1.  INTRODUCTION 

In  his  review  of  Ambient  Noise  in  the  Sea,  Urick  (1984)  commented  on 
the  prolific  nature  of  literature  concerning  theories  of  sound 
generation  at  the  sea  surface  and  the  measurements  of  the  temporal  and 
spatial  spectral  characteristics.  The  idealized  spectra  suggested  by 
Urick  was  in  close  agreement  with  the  schematic  proposal  by  Wenz 
(1962).  The  spectral  characteristics  for  frequencies  greater  th?.i  500 
Hz  are  consistent  with  the  observations  of  Knudsen  (1948)  and  V.enz's 
"rule  of  fives"  that  the  noise  level  can  be  described  by 

NL(f,U)  =  25  -  5L0G(  f )  +  5LOG(U/5)  [dB  re  20  v’.Vm2 )  , 


*  Kildare  Corporation,  95  Trumbull  Street,  New  London,  CT  06320  USA. 


where  f  is  frequency  (kHz)  and  U  is  wind  speed  (knots).  However,  as 
Wenz  observed  in  the  10-  to  500-Hz  band,  the  measured  noise  levels 
were  often  variable  and  dominated  by  shipping  noise.  The  shape  of  the 
spectrum  was  also  found  to  vary  from  a  positive  slope  to  a  steep 
negative  slope. 

Kerman  (1984)  showed  that  "the  amalgamated  observations  of  the 
ambient  noise  reveal  a  similarity  structure,  both  in  the  acoustical 
spectrum  and  wind  dependency."  For  frequencies  greater  than  the 
local  maximum  in  the  300-  to  500-Hz  range,  Kerman  found  that  the 
normalized  measured  spectral  characteristic  was  proportional  to 
f“2(6  dB/octave).  Furthermore ,  he  showed  that  the  noise  intensity 
was  proportional  to  the  cube  of  the  friction  velocity  (u^)  prior  to  a 
critical  friction  velocity  (u^)  which  is  determined  by  the  minimum 
phase  velocity  of  the  gravity-capillary  waves.  Wave  breaking  was 
associated  with  this  critical  condition,  and  for  u^  >  ux.  the  noise 
intensity  was  found  to  increase  with  u£-->.  These  observations  were 
found  to  be  consistent  with  a  large  number  of  experimental  observations 
cited  by  Kerman.  He  speculated  that  the  two  observed  regions  of  ambient 
noise  wind  speed  dependency  indicated  the  presence  of  two  sound  source 
generation  mechanisms  or  one  mechanism  that  changes  sensitivity.  Never¬ 
theless,  since  breaking  waves  are  known  to  produce  bubbles,  spray, 
splash,  and  turbulence,  combinations  of  these  mechanisms  may  explain 
the  production  of  sound  at  frequencies  >500  Hz.  Kerman  also  observed  a 
variability  in  the  region  <500  Hz. 

At  frequencies  <500  Hz,  it  has  been  difficult  to  acquire  ambient 
noise  data  that  are  attributable  solely  to  local  effects.  The  reason 
for  this  is  the  dominance  of  shipping  noise  in  the  Northern  Hemisphere, 
the  problems  of  self-noise,  the  size  required  for  directional  measure¬ 
ment  systems,  and  the  relatively  low  long  range  propagation  loss. 

Reviews  of  low-  to  mid-frequency  (20  to  1000  Hz)  ambient  noise 
measurements  (Carey  (1985,1988))  not  dominated  by  shipping  at 
frequencies  <500  Hz  showed  the  presence  of  locally  wind-generated 
noise.  Wittenbom  (1976)  showed  a  wind  speed  dependence  with  a  linear 
velocity  dependence  prior  to  u~6  m/s  and  a  quadratic  dependence  at 
greater  wind  speeds.  His  results  were  not  based  on  actual  local 
measurements  of  environmental  variables  such  as  wind  speed.  Measure¬ 
ments  performed  with  vertical  arrays  in  the  sparsely  shipped  Southern 
Hemisphere's  Fiji  Basin  (Bannister  (1981),  Burgess  and  Kewley  (1983), 
Browning  (1986))  were  also  found  to  have  low  frequency  wind  noise  with 
two  distinct  wind  speed  dependencies.  When  compared  to  other  shallow 
and  deep  water  noise  measurements,  these  results  led  to  the  conclusion 
that  at  the  lower  frequencies  sound  was  generated  near  the  sea  surface 
with  one  wind  speed  dependent  mechanism  prior  to  wave  breaking  and 
another  after  wave  breaking  (Carey  (1985,1988)). 

Two  interesting  experiments  have  recently  been  performed  featuring 
visual  as  well  as  acoustic  identification  of  wave  breaking  events. 
Hollett  (1988)  performed  an  experiment  in  the  Mediterranean  Sea  with  a 
vertical  array  of  hydrophones  (3  nested  apertures  (32  phones  each)  with 
center  frequencies  at  375,  750,  and  1500  Hz)  and  simultaneous  video 
recording  of  the  sea  surface  area  intersected  by  the  endfire  beam  of 
the  array.  Figure  1  shows  the  spectral  events  which  result  from  a 
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large  breaker.  The  spectra  shown  have  not  been  corrected  for  the 
prewhitening  of  the  data;  i.e.,  the  lower  frequency  spectral  content  is 
more  pronounced  than  shown.  The  breaking  event  occurs  between  3  to  7 
seconds  and  is  seen  to  be  a  random  collection  of  spectral  peaks  most 
pronounced  in  the  low  frequency  (<300  Hz)  range. 

Farmer  (1989)  performed  an  experiment  in  200  m  of  water  with  a 
hydrophone  14  m  below  the  surface  and  was  able  to  identify  the 
occurrence  of  wave  breaking  by  examining  the  video  obtained  from  a 
subsurface  camera.  The  simultaneously  recorded  acoustic  data  were 
examined  and  found  to  show  that  the  breaking  waves  radiate  sound  to 
frequencies  as  low  as  50  Hz. 


Figure  1.  Spectra  of  a  wave  breaking  as  measured  by  Hollett  (1988)  are 
shown  as  a  function  of  time.  The  lower  frequency  spectral  estimates 
have  not  been  corrected  for  pre-emphasis.  This  collection  of  acoustic 
spectra  shows  the  wave  breaking  event  to  be  a  random  collection  of 
spectral  peaks. 

The  Hollett  and  Farmer  results  clearly  show  that  breaking  waves 
are  a  source  of  low  frequency  sound  and  that  at  frequencies  <  500  Hz 
the  ambient  noise  spectrum  in  the  absence  of  shipping  is  wind  speed 
dependent.  These  results  are  important  with  respect  to  measurements 
performed  with  vertical  and  horizontal  arrays  where  the  directional 
noise  properties,  in  the  absence  of  shipping,  are  due  to  the 
interaction  of  wind-generated  noise  and  the  basin  boundaries. 


Since  breaking  waves  produce  sound  at  frequencies  <  500  Hz,  a 
question  concerning  the  generating  mechanisms  naturally  arises.  This 
paper  shall  examine  the  role  of  spray,  impact,  and  the  collective 
oscillations  of  bubble  plumes  and  clouds  proposed  by  this  author  and 
Prosperetti  (1985).  First,  the  properties  of  low  frequency  ambient 
noise  in  terms  of  the  wind  speed  dependence  and  spectral  characteristic 
will  be  presented.  This  section  will  feature  recent  experimental  work 
since  the  Lerici  meeting  (Carey  (1988)).  Second,  the  properties  of 
mixtures  will  be  discussed  and  theoretical  expressions  for  spherical 
and  cylindrical  bubble  clouds  presented.  Calculations  will  be  presented 
for  the  noise  production  in  a  hypothetical  wave  breaking  event  based  on 
an  evolutl-n  modej.  due  to  Monahan  that  demonstrates  that  bubble  plume 
and  cloud  oscillation  are  candidate  mechanisms  for  the  production  of 
noise. 


2.  A  REVIEW  OF  RECENT  LOW  FREQUENCY  NOISE  MEASUREMENTS 

Since  a  voluminous  amount  of  ambient  noise  literature  exists  (Urick 
(198A),  Kerman  (1984),  (1988)),  the  purpose  of  this  review  is  to  simply 
concentrate  on  ambient  noise  characteristics  <500  Hz.  This  characteri¬ 
zation  will  be  based  on  the  determination  of  two  parameters,  the 
spectral  slope  m(f)  and  the  wind  speed  dependency  n(f)  defined  by 

ANL(f ,f0,U,U0)  =  20n(f )L0G(U/Uo)  +  10m(f ) *LOG(f /f0) , 

where  ANL  is  the  difference  in  measured  noise  level  at  frequencies  f 
and  f0  and  wind  speeds  U  and  UQ.  These  parameters  are  usually  estimated 
by  the  use  of  regression  analysis  and,  in  some  instances,  by  a  visual 
least-squares  estimation.  We  have  previously  stated  that  the  measure¬ 
ment  of  low  frequency  wind-dependent  noise  is  difficult  due  to  the 
presence  of  shipping  noise,  self-noise,  and  sound  propagation 
conditions.  For  this  reason,  we  have  been  selective  with  respect  to 
those  measurements  included  in  this  summary. 

2.1  The  Low  Frequency  Spectral  Slope 

In  a  previous  paper  (Carey  (1988)),  the  experimental  results  from 
Wittenbom  (1976)  were  presented  as  an  example  of  deep  ocean  wind¬ 
generated  ambient  noise.  This  experiment  consisted  of  a  vertical 
string  of  omnidirectional  hydrophones,  one  of  which  was  positioned 
below  the  sound  channel’s  critical  depth  to  minimize  the  hydrophone's 
reception  of  distant  generated  noise.  Additional  results  from  this 
experiment  may  be  found  in  a  recent  article  by  Shooter,  De  Mary,  and 
Wittenbom  (1990).  The  hydrophone  data  were  remotely  recorded  on  a 
time-indexed  magnetic  tape.  Wind  speed  estimates  were  based  on 
interpolating  weather  reports  of  ships  transiting  the  area.  Stationary 
periods  were  selected  by  examining  the  continuous  recording  of  the  300- 
to  500-Hz  band  of  the  deep  hydrophone.  However,  the  lack  of  concurrent 
wind  speed  and  air-sea  temperature  difference  is  a  major  shortcoming  of 
this  experiment.  Nevertheless,  the  results  clearly  show  wind-dependent 
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noise  characteristics  between  20  and  500  Hz.  The  15-knot  curve  shown 
in  Fig.  2  has  spectral  slopes  of  m=0, 1,1/3.  The  corresponding  wind 
speed  factor  (n)  determined  by  the  difference  in  noise  levels  at  12  and 
15  knots  and  a  frequency  of  100  Hz  was  n=1.65. 


Figure  2.  Source  spectrum  level  (SL)  and  noise  spectrum  level  (NL)  are 
shown  on  this  plot  of  spectrum  level  vs  frequency  for  wind  speeds  >13 
knots  for  several  diverse  oceanographic  areas. 

Several  investigators  have  used  vertical  arrays  to  measure  the 
locally  generated  noise  in  deep  water.  Most  notable  are  the  works  of 
Burgess  and  Kewley  (1983)  in  the  South  Fiji  Basin  and  Australian  waters 
and  Kennedy  (1990)  in  the  Tongue  of  the  Ocean  in  the  Bahamas.  In 
addition  to  the  use  of  steerable  arrays,  these  experiments  were 
conducted  in  areas  to  minimize  shipping  noise  and  had  concurrent 
measurement  of  wind  speed. 

The  Burgess  and  Kewley  (1983)  experiment  used  a  180-m  steerable 
array  at  a  depth  of  300  m  in  deep  water  to  avoid  mixed  layer  ducting 
effects.  Both  upward  and  downward  looking  beams  were  used  to  determine 
the  sea  surface  source  level  (SL) ,  dB  re  1  yPa2/(sr«Hz) .  Since  their 
original  publication,  these  authors  have  re-examined  their  linear  re¬ 
gression  analysis  of  the  wind  speed  dependency  and  have  concluded  that 
two  wind  speed  dependent  regions  exist,  i.e.,  one  prior  to  and  one 
after  the  onset  of  wave  breaking  with  n~1.5  (Kewley  et  al.  (1990)). 
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Their  frequency-dependent  source  level  estimate  is  shown  in  Fig.  2.  In 
examining  Fig.  2,  one  needs  to  differentiate  between  noise  level  (NO 
and  source  level  (SL)  estimates.  We  have  chosen  to  plot  these  results 
on  the  same  scale  to  show  that  there  is  no  strong  spectral  dependence 
at  wind  speeds  greater  than  6  m/s  (12  knots).  In  the  case  of  the  Kewley 
(1990),  we  observe  at  the  most  m  =  +1/3  between  100-200  Hz. 

The  Kennedy  (1990)  experiment  was  performed  with  a  measurement 
array  of  seven  octavely  nested  four  wavelength  linear  apertures 
covering  the  40-  to  4000-Hz  band.  Meteorological  measurements  were 
performed  on  Andros  Island.  In  his  experiment,  Kennedy  was  able  to 
model  the  propagation  in  the  Tongue  of  the  Ocean.  Since  we  know 
(Talhara  (1964))  that  the  vertical  distribution  of  ambient  noise  N(6) 
can  be  written  as 

N(0)  =  dl/dQ  =  DPWg(6' )exp(-2ar)/cos0' (l-By  exp(-4ar)), 

where  W  =  tsinO'/sinO]  =  [Cs/Cr]2, 

D  =  source  density  at  surface  (number/m^), 

P  =  power  (watts/sr), 

dl/dfl  =  intensity/unit  solid  angle  (watts/m^/sr) , 
a  =  frequency-dependent  absorption  coefficient, 

Y  =  surface  reflection  coefficient, 

B  =  reflection  coefficient  of  the  bottom, 

0*  =  source  angle, 

0  =  angle  at  receiver,  and 

Cg,Cr  =  speed  of  sound  at  source  and  receiver, 

we  see  immediately  that,  given  a  knowledge  of  the  environmental  factors, 
one  could  estimate  g(0')  given  a  measurement  of  N(0) .  Kennedy 
measures  N(f,0)  and  determines  with  a  similar  but  more  complicated 
method  of  curve  fitting  the  best  fit  to  D»P*g(0’).  Kennedy's 
results  are  shown  in  Fig.  2  for  a  wind  speed  of  8  m/s  (16  knots).  These 
results  show  a  practically  white  spectral  curve.  The  novel  aspect  of 
Kennedy's  experiment  is  his  ability  to  estimate  the  source  directional 
characteristic  g(0')  as  a  function  of  wind  speed  and  frequency.  Shown 
in  Fig.  3  are  two  examples  of  N(f,0)  from  this  estimation  based  on 
his  measured  data.  Figure  3A  shows  the  estimated  pattern  at  low  wind 
speeds  prior  to  whitecaps  being  present.  A  broad  maximum  is  observed  at 
600  Hz  that  is  directional.  At  the  lower  frequencies,  less  direction¬ 
ality  is  observed.  Kennedy  found  that  with  whitecaps  present  (Fig.  3B) 
the  vertical  directional  spectra  were  consistent  with  a  surface 
distribution  of  point  dipole  sources  over  the  entire  frequency  range 
with  a  broad  maximum  at  400  Hz . 

The  structure  shown  in  Fig.  3A  can  also  be  modeled  by  point 
dipoles  near  the  broad  maximum  at  600  Hz,  but  at  the  lower  frequencies 
(<300  Hz)  the  structure  is  consistent  with  volume-distributed  monopole 
sources  or  noise  from  a  distance.  Thus,  Kennedy’s  results  not  only  show 
a  slowly  varying  spectrum  level  with  frequency  but  also  a  dipole 
characteristic  for  sound  generated  near  the  sea  surface  associated  with 
wave  breaking. 
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Figure  3.  Estimated  noise  spectrum  level  vs  elevation  angle  as 
determined  by  Kennedy  (1990)  from  comparison  of  the  field  from  surface 
distributed  dipoles  and  measured  results. 


Marrett  and  Chapman  (1989)  have  performed  measurements  of  low 
frequency  (15  to  250  Hz)  ambient  noise  in  the  South  Fiji  Basin  with  a 
towed  line  array.  Local  wind  speed  measurements  were  performed.  Shown 
in  Fig.  2  are  the  25-knot  noise  level  results.  The  spectrum  again  has 
a  fairly  white  character.  For  wind  speeds  greater  than  15  knots,  the 
wind  speed  dependency  factor  was  found  to  be  n=1.32  at  250  Hz  to 
n*2.13  at  30  Hz. 

These  recent  deep  water  ambient  noise  measurements  show  that  for 
wind  speeds  greater  than  6  m/s,  the  sea  surface  sound  has  a  slowly 
varying  spectral  shape  with  a  broad  maximum  ranging  from  300  to  500  Hz. 
In  addition,  we  are  shown  in  Fig.  2  the  results  from  three  shallow 
water  experiments.  Ferla  and  Kuperman  (1984)  have  reported  results  on 
the  wind  generated  source  spectrum  levels  for  wind  speeds  between  10  to 
40  knots.  Ferla  and  Kuperman  used  propagation  loss  measurements  and 
environmental  information  in  a  wave  theoretic  noise  model  to  derive 
these  source  levels  from  measured  noise  levels.  As  can  be  seen,  their 
results  are  consistent  with  the  deep  water  curves.  Also  shown  are 
results  from  Wille  (1984)  in  the  North  Sea  at  40  to  50  knots  and 
Hollingberger  in  an  Alaskan  fjord  at  14-15  knots.  To  obtain  the  curve 
attributed  to  Wille,  we  have  applied  an  estimated  but  relative 
correction  factor  based  on  measured  propagation  loss  to  their  measured 
noise  levels.  The  results  art  all  similar  in  shape  although  different 
in  level. 

This  characteristic,  the  broad  maximum  between  300  to  500  Hz,  has 
been  recognized  in  wind-generated  ambient  noise  for  quite  some  time 
(Wenz  (1962),  Piggott  (1964)),  but  not  identified  with  such  a  diverse 
set  of  experiments.  The  broad  maximum  that  is  seen  between  300  to  500 
Hz  in  most  of  the  curves  shown  in  Fig.  2  shows  a  decreasing  spectrum 
level  as  frequency  decreases.  This  decrease  corresponds  to  a  spectral 
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slope  factor  of  nr-1/3  or  1  dB/octave.  This  characteristic  has  often 
been  compared  to  the  relative  spectrum  of  noise  due  to  a  spray  of  water 
droplets  developed  by  Franz  (1959).  However  the  Franz  spectrum  has  a 
steeper  (1.3-1. 4  dB/active  octave)  slope  than  those  shown  here;  never¬ 
theless,  the  observation  of  their  similarity  is  worthy  of  detailed 
consideration. 


2.2  The  Low-Frequency  Wind  Speed  Dependency 


In  the  previous  discussion  of  low-frequency  noise  measurements,  we  have 
continually  indicated  that  the  wind  speed  dependency  characterized  by 
the  factor  n(f)  has  been  observed  to  range  from  1.3  to  2.5.  Wind  speed 
has  traditionally  been  used  to  provide  an  index  to  the  level  of  noise 
to  be  expected  at  the  higher  frequencies.  But  even  at  high  frequencies, 
the  wind  speed  dependency  factor  may  be  variable  and  complex.  To 
illustrate  the  complexity  of  wind  speed  dependency  factor  n,  we  have 
produced  the  schematic  shown  in  Fig.  4,  which  has  been  adapted  from  the 
work  of  Wille  and  Geyer  (1984).  These  data  were  obtained  in  the  North 
Sea,  a  region  with  appreciable  shipping-generated  noise.  The  spec  .'ic 
data  set  illustrates  what  we  believe  is  found  in  the  licerature 
concerning  the  wind  speed  dependency.  The  spectral  curve  at  20  kHz 
below  a  wind  speed  class  of  10  shows  no  dependency  (n=0);  between  wind 
speed  classes  11  and  13,  a  dependency  of  n=2;  and  at  higher  wind 
speeds,  a  reversal  of  slope  and  a  decrease  in  noise  level. 


WIND  CLASS 


WIND  WIND 

CLASS  SPEED 
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Figure  4.  Spectrum  level  of  ambient  noise  vs  wind  speed  class  in  the 
North  Sea  as  measured  by  Wille  (1984).  These  spectra  show  the  general 
wind  speed  dependency  of  ambient  noise. 


This  trend  is  seen  to  persist  in  some  degree  down  to  5  kHz.  The 
decrease  in  sound  level  at  wind  speeds  greater  than  10  m/s  has  been 
attributed  to  bubble  absorption  and  sound  trapping  in  the  near-surface 
sound  duct  produced  by  sound  speed  gradients  resulting  from  the  presence 
of  microbubbles  (Wille  and  Geyer  (1984),  Farmer  and  Vagle  (1989)).  At 
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1  kHz  we  observe  a  region  of  no  wind  speed  dependence  below  a  wind  speed 
class  of  12  and  a  wind  speed  dependency  of  n=1.75  at  greater  wind 
speeds.  At  lower  frequencies,  =200  Hz,  a  wind  speed  dependency  is 
only  observed  at  the  higher  wind  classes  (W.C.>16,  W.S.>12.5  m/s). 

In  general  what  is  observed  in  measurements  of  wind  speed 
dependence  is  a  noise-limited  region,  a  transition  region,  and  a  high 
wind  speed  region.  These  regions  are  shown  in  Fig.  5  for  representative 
wind  speeds  and  values  of  n.  Classifying  measurements  i^  these  regions 
may  aid  in  the  interpretation  of  low  frequency  noise  levels  and  explain 
the  variation  found  in  the  published  values  of  n(f). 


WIND  SPEED  (m/s) 

Figure  5.  Wind  speed  dependency  (n(f))  vs  wind  speed  for  frequency 
<500  Hz.  Three  regions  are  identified,  a  noise  limited  region,  a 
transition  region,  and  a  high  wind  speed  region. 

We  have  used  this  scheme  to  examine  the  value  for  n(f)  for  noise 
results  obtained  at  wind  speeds  greater  than  6.5  m/s  on  selected 
experimental  results.  The  summary  of  this  analysis  is  shown  in  Fig.  6 
along  with  the  values  reported  by  Crouch  and  Burt  (1972),  Piggott 
(1964),  and  Marrett  and  Chapman  (1990).  It  is  important  to  mention  that 
all  results  were  used  in  our  summary  curve  shown  in  Fig.  6.  The  trend 
is  clear  for  high  winds  speeds;  n(f)  has  a  constant  value  ~2.5  until 
a  frequency  of  50  Hz  and  ~1 . 2  value  for  frequencies  greater  than  300 
Hz.  The  variation  of  n(f)  with  f  between  50  and  200  Hz  appears  to  be 
real.  Kerman  (1984)  has  observed  that  all  high  frequency  results  (>500 
Hz)  yield  a  consistent  set  of  characteristics.  This  indicates  that 
sound  production  at  these  higher  frequencies  may  be  from  the  same 
mechanisms,  such  as  bubble,  splash,  and  spray.  However,  the  variation 
of  n(f)  at  the  lower  frequencies  may  indicate  a  different  mechanism. 

2.3  Low-Frequency  Noise  Characteristics 

We  now  have  an  answer  to  the  question  "What  are  the  characteristics  of 
low  frequency  (<500  Hz)  ambient  noise?"  First,  we  cited  evidence  that 
showed  breaking  waves  produce  sound  with  frequency  content  as  low  as 
30  Hz.  The  Hollett  spectra  appear  to  be  a  random  collection  of  spectral 
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Figure  6.  Wind  speed  dependency  factor  n(f)  vs  frequency  for  wind 
speeds  >13  knots.  Shown  are  the  results  of  Crouch  and  Burt  (1972), 
Piggott  (1964),  and  Marrett  and  Chapman  (1990)  and  an  average  of  all 
data . 

peaks  spread  in  frequency  and  time  during  the  breaking  event.  These 
results  are  similar  to  those  of  Farmer  (1989),  although  in  some 
instances  Farmer  observed  a  broadband  event  representative  of  impact 
noise  at  the  beginning  of  the  wave  breaking,  followed  by  the  subsequent 
random  collection  of  spectral  peaks. 

Second  we  showed  that  measurements  of  local  wind-generated  noise 
at  wind  speeds  greater  than  6  m/s  have  a  broad  maximum  between  300  and 
500  Hz.  The  position  of  this  maximum  shifts  to  lower  frequencies  as  the 
wind  speed  increases.  The  decrease  in  spectrum  level  with  a  decrease  in 
frequency  below  this  broad  maximum  is  less  than  m=l/3  (1  dB/octave) . 

Third,  we  showed  that  the  wind  speed  dependency  (n(f))  of  local 
wind  generated  noise  is  a  complex  function  and  may  be  characterized  by 
three  regions.  In  the  high  wind  speed  case,  the  low-frequency  noise 
level  (dB  re  1  yPa^/Hz)  increased  with  a  wind  speed  dependency 
n(f)=1.5  at  200  Hz,  which  ranged  between  1.3  at  500  Hz  to  2.5  at  30  Hz. 

Finally,  the  results  of  Kennedy  (1990)  clearly  show  that  low 
frequency  sound  produced  at  the  sea  surface  after  the  occurrence  of 
wave  breaking  has  a  dipole  characteristic.  Thus,  the  problem  is 
reduced  to  finding  a  sound-generating  mechanism  that  produces  a 
broadband  event  followed  by  a  random  collection  of  spectral  events, 
which  scale  according  to  power. 

2.4  The  Whitecap  Index 

Wilson  (1983)  proposed  that  above  200  Hz,  the  noise  levels  vary  in 
proportion  to  the  "whitecap"  index,  W(u),  of  Ross  and  Cardone  (1979). 
Wilson  proposed  three  regions  of  wind  speed  dependency: 


10  100  1000 
FREQUENCY  (Hz) 
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I . 

W(u)  =  0, 

u  <  4.5  m/s 

II . 

W(u)  =  (4 . 6xlO-3 )U3  -  (4.9x10' 

2)U2 

4-  (4 . 63xl0_1  )U  -  1.5, 

4.5  <  U  <  15  m/s 

Ill . 

W(u)  =  (20.97MU/15)1-5, 

15  m/s  <  U  . 

The  first  region  has  no  whitecaps;  the  second  was  attributed  to  an 
empirical  fit  to  the  data  of  Ross  and  Cardone;  the  third,  or  high  wind 
speed  region,  was  based  on  Wilson's  analysis  of  Perrone's  data.  The 
speculation  put  forth  by  Wilson  was  the  relationship  between  whitecaps 
and  acoustic  noise  was  in  three  individual  wind  speed  dependent  regions. 

We  have  observed  a  wind  speed  dependency  of  n(200)~1.5, 
a  noise  intensity  variation  of  However,  if  low-frequency  sound  is 

associated  with  the  breaking  waves  characterized  by  the  whitecap  index, 
which  vary  as  ,  then  why  does  the  low-frequency  noise  level  only 
vary  with  and  why  is  this  factor  so  variable.  Kerman  (1984)  has 

developed  an  explanation  for  the  higher  frequency  wind  speed 
dependence  based  on  bubbles  with  radii  comparable  to  the  Kolmogorov 
length  scale,  but  this  explanation  may  not  apply  at  these  frequencies. 

Wu  (1980)  stated  that  the  whitecap  index  W(u)  should  be  related 
to  the  energy  flux  of  the  wind  under  equilibrium  conditions.  The 
energy  flux  (E)  or  the  rate  of  doing  work  is  related  to  the  wind 
stress  (t)  and  a  surface  drift  current  (V), 

W(U10)  «  E  =  tV  <*  tU*  <x  2  u3q  , 

where  W(U^q)  is  the  percentage  of  the  sea  surface  covered  by  white 

caps,  t  is  the  shear  stress  at  the  surface,  U*is  the  friction  velocity, 
and  U10  is  the  10-meter-elevation  wind  speed.  This  relationship 
between  the  white  cap  index  and  the  wind  stress  coefficient  (C^q)  ties 
the  whitecap  index  to  fundamental  parameters  governing  the  exchange  of 
momentum,  mass  and  energy  in  the  sea  surface  interaction  zone. 

Kitaigorodskii  (1972)  outlined  the  fundamental  dependency  of  the 
wind  stress  coefficient  on  three  quantities 

C(Z,U)  -  C(Z/hs,Ri,Res)  , 

where  Z  is  the  observation  point,  hg  is  the  characteristic  scale  of 
the  roughness,  R^  is  the  Richardson  number,  and  Res  is  the  Reynolds 
number  for  roughness.  Kitaigorodskii  shows,  for  a  moving  boundary  under 
neutral  stability  conditions,  that 

U(Z)/UA  =  1 / k  ln(Z/hs)  =  C(Z,U)-1/2, 

where  U*=  \/Tt/p)  is  the  friction  velocity,  k  is  Von  Ka man's  universal 
constant  and  hs  is  given  by 

®  1/2 

hs  -  (2  /  S(<j)exp(- 2Kg/ojUA)dw) 
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with  S(u)  being  the  frequency  spectrum  of  the  wave  field.  In 
general,  the  wind  stress  coefficient  is  dependent  on  the  stability  of 
the  boundary  layer  as  well  as  the  roughness  scale  modified  by  the 
motion.  Given  that  the  whitecap  index  is  proportional  to  this  quantity 
to  the  3/2  power,  we  can  expect  variability. 

The  whitecap  index  can  be  approached  by  either  measuring  the 
index  directly  or  by  measuring  the  wind  stress  coefficient.  Monahan 
(1990)  (see  also  O'Murcheartaigh  and  Monahan  (1986),  Monahan  and 
O'Murcheartaigh  (1981))  has  determined  the  wind  speed  variation  of  the 
whitecap  index  by  fitting 

X 

W(Uio)  =  aU1Q 

to  measured  fractional  whitecap  coverage  data.  Monahan  finds  that 
X  =  3.41  provides  the  best  fit  to  all  data  sets,  but  individual  sets  of 
data  yielded  values  of  X  between  2.55  and  3.75.  Monahan  has  also 
proposed  classifying  whitecaps  as  class  A,  young,  and  class  B,  old. 
Wu(1981 , 1986 )  contends  that  C(U10)  «  an(j  consequently  a  variation 

of  W  with  1)3-75.  However,  Wu  (1980)  in  a  revisit  to  this  question 
recognized  the  linear  dependence  of  C(U^q)  on  wind  speed.  His  review  of 
measurement  of  the  wind  stress  coefficient  as  a  function  of  U^m  showed 
most  measured  values  of  m  between  a  value  of  1  and  1.3.  (Also,  see  Large 
(1981),  Donelan  (1982)  and  Smith  (1980).)  The  reader  is  cautioned 
concerning  the  use  of  relationships  between  the  whitecap  index  and 
C(U10). 

Amorocho  and  De  Vries  (1980)  examined  the  issue  of  variability 
by  plotting  C(U10)  and  u*  versus  U^q  for  a  wide  range  of  wind  speed 
conditions.  The  results  of  their  curve  fits  are  shown  in  Fig.  7.  In 
general,  they  find  three  distinct  wind  speed  regions.  The  first  region 
is  found  prior  to  the  onset  of  breakers  with  C^o  constant.  The  second 
region,  labeled  as  a  transition  region,  is  for  wind  speeds  between  7 
and  20  m/s.  In  this  region,  both  U*  and  C^q  have  linear  dependencies 
on  wind  speed.  This  region  would  correspond  to  W(u)  «  u^50(  which  is 
larger  than  Monahan's  observations.  Finally,  for  wind  speeds  greater 
than  20  m/s,  a  saturation  region  with  C^q  again  constant  and  an 
expected  variation  of  W(U^q)  * 

These  results  are  remarkably  similar  to  the  acoustic  noise  level 
characteristics.  Since  most  observations  of  noise  are  in  the  7-  to 
20-m/s  wind  range,  we  can  expect  a  variation  in  our  wind  speed 
dependency  similar  to  the  variation  in  wind  stress  coefficient.  Since 
we  also  expect  a  dependency  of  the  wind  stress  coefficient  on  the 
Richardson  number  and  the  sea  state  spectrum,  ambient  noise  levels  may 
also  be  affected.  Thus,  one  may  conclude  that  ambient  noise  measure¬ 
ments  should  be  performed  that  concurrently  measure  the  wind  speed 
(U10) ,  the  air-water  temperature  difference  as  a  simplified 
measure,  the  sea  state  spectrum  as  an  indicator  of  roughness,  and  the 
moisture  content  necessary  for  correct  estimation  of  Re. 
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Figure  7.  The  friction  velocity,  U*,  and  wind  stress  coefficient,  C^q» 
vs  wind  speed  taken  from  Amorocho  and  DeVries  (1980).  These  regions  are 
identified  as  a  constant  slope  region  prior  to  wave  breaking,  a 
transition  region,  and  a  constant  slope  after  breaker  saturation. 


3.  BREAKING  WAVES:  BUBBLE  PLUMES  AND  CLOUDS 

The  association  of  wave  breaking  with  the  production  of  low  frequency 
sound  and  subsequently  deep  ocean  ambient  noise  characteristics  leads 
us  to  examine  bubble  production  from  this  wave  breaking  process.  Fresh 
water  laboratory  measurements,  such  as  the  experiments  by  Ming-Yang  Su 
(1984)  with  three-dimensional  waves,  show  a  typical  sequence  of  events 
from  the  initial  wave  breaking  to  the  formulation  of  a  columnar  bubble 
plume  to  a  more  diffuse  cloud  swept  by  the  subsurface  orbital  motion 
and  finally  to  dispersion  of  near-surface  bubbly  features.  While  fresh 
water  experiments  may  be  useful  in  visualizing  the  sequence,  the  bubbly 
mixture  in  these  types  of  experiments  appears  to  be  composed  of  larger 
bubbles,  faster  rise  times,  and  consequently  different  acoustic  char¬ 
acteristics.  The  difference  in  the  bubble  sizes  and  rise  times  may  also 
result  in  different  sound  production  mechanisms. 

Thorpe  (1986,1982)  and  his  colleagues  have  performed  a  series  of 
experiments  with  upward- looking  sonars  in  lakes  and  in  the  deep  ocean. 
He  observed  bubble  clouds  convected  to  meter  depths  initially  by  wave 
breaking  turbulence  and  vorticity  and  subsequently  ordered  patterns 
consistent  with  Langmuir  circulation.  Thorpe's  results  show  an 
exponential  decrease  in  volume  scatter  strength  with  depth,  a  mean 
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depth  of  penetration  that  was  proportional  to  U^q  and  the  air-water 
temperature  difference  (A0^w)  and  finally  a  different  cloud 
characteristic  for  stable  (A0^W>O)  and  unstable  (A0^W<O)  conditions. 

For  example,  the  stable  condition  was  found  to  result  in  a  "billowy” 
cloud  structure,  whereas  the  unstable  condition  resulted  in  a  columnar 
characteristic.  Recently,  observations  by  Crawford  and  Farmer  (1987) 
confirmed  these  effects  as  well  as  an  exponential  distribution  of  bubble 
density  with  an  "e  folding"  depth  between  0.7  and  1.5  m,  a  near-surface 
bubble  density  variation  with  u^o^-'^>  a  weak  dependence  on  the 
Langmuir  circulation,  and  "v”-shaped  columnar  clouds  for  (A0^W<O). 

Thorpe  also  observed  pronounced  differences  in  clouds  produced  by 
wave  breaking  in  fresh  and  salt  water.  Thorpe  attributed  these  differ¬ 
ences  to  chemical  effects  discussed  by  Scott  (1975)  to  explain  Monahan’s 
(1969,  1971)  observation  concerning  fresh  water  white  caps.  That  is, 
under  nearly  identical  physical  conditions,  bubble  distributions  pro¬ 
duced  in  salt  water  have  a  smaller  mean  radii  and  a  larger  number  of 
bubbles.  According  to  Scott  (1986),  surface  chemical  effects  can  be  an 
important  factor  preventing  coalescence,  and  "Significant  differences 
observed  in  fresh  water  and  salt  water  white  caps  may  be  ascribed  to 
these  effects,  bubbles  in  salt  water  being  (as  a  result)  greater  in 
number,  smaller,  more  densely  packed,  carried  deeper,  and  slower  to 
rise  to  the  surface  than  those  formed  in  fresh  water  by  a  similar  wave 
breaking  event."  Pounder  (1986)  also  showed  a  distinct  temperature- 
dependent  difference  between  distilled  water  (coalescence  occurs)  and 
salt  water  (coalescence  does  not  occur).  He  attributed  this  difference 
to  ionic  effect.  However,  Pounder's  laboratory  observations  support  the 
conclusion  drawn  by  Scott  and  Thorpe  that  microbubble  distributions 
result  from  salt  water  wave  breaking. 

Bubble  size  measurements  (made  by  Medwin  (1977),  Kalovayev  (1976), 
Johnson  and  Cooke  (1979),  Bouguel  (1985),  and  Crawford  (1987)  and  then 
reviewed  by  McDaniel  (1987))  yield  distributions  with  most  probable 
bubble  radii  between  50-70  pm,  with  an  exponential  decrease  in  number 
with  an  increase  in  radii.  Kalovayev  observed  that  at  wind  speeds  of 
13  m/s,  all  bubbles  were  less  than  350  ym,  although  the  position  of 
the  distribution  maxima  was  observed  to  shift  to  larger  radii  with 
increasing  wind  speed  and  depth. 

Recently,  Monahan  (1990,1987)  has  developed  a  hypothetical  model 
for  the  evolution  of  a  bubble  plume  and  cloud  from  a  breaking  wave 
consistent  with  the  results  discussed  here.  Shown  in  Fig.  8  are  the 
results  of  his  analysis  for  a  wind  speed  of  13  m/s.  The  second  panel 
shows  his  "8"  plume,  a  bubble  size  distribution  derived  from  his 
aerosol  generation  model,  and  a  bubble  size  distribution  based  on  the 
measurements  of  Johnson  and  Cooke  (1979). 

Monahan  uses  exponential  variation  in  depth,  cross  sectional  area, 
and  time  with  measured  "e  folding"  characteristics  to  scale  the  results 
for  these  a,  8,  and  y  plumes.  The  a  plume  occurs  within  1-2  s  of  the 
breaking  event  and  its  characteristic  depth  of  0.5  m  and  volume  fraction 
(v.f . )  of  4-8x10“2  are  based  on  the  extrapolated  bubble  size  distribution 
shown  for  the  8  plume.  The  8  plume  is  estimated  to  have  a  duration  between  1 
and  10  s  and  a  v.f.  based  on  an  integrated  size  distribution  of  l-2xl0~4 . 
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Third,  che  characteristic  of  the  cloud  has  a  bubble  size  distribution 
consistent  with  the  measurements  of  Johnson  and  Cooke  and  a  scaling  based  on 
cloud  sea  surface  ocean  to  white  cap  area  of  25:1.  The  duration  of  this 
feature  is  estimated  to  be  on  the  order  of  100  s  and  an  average  v.f.  between 
10-^  to  lO-^.  Monahan’s  estimates,  although  carefully  derived,  are  still 
very  uncertain  and  illustrate  that  the  measurement  of  v.f.  is  required  to 
characterize  the  wave  breaking  process. 
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R,  BUBBLE  RADIUS,  IN  Mm 


Figure  8.  Evolutionary  model  of  a  wave  breaking  bubble  plume  and  cloud 
shown  as  a  sequence  of  snapshots  covering  the  time  periods  of  1  s,  10  s, 
100  s.  Estimates  of  v.f.  and  v.f.  distributions  are  from  Monahan  (1988). 


The  measurement  of  the  bubble  size  distribution,  void  fraction,  and 
spatial  characteristics  of  bubble  clouds  and  plumes  has  been  shown  to  be 
difficult.  Nevertheless,  general  characteristics  are  beginning  to  emerge 
concerning  the  production  of  bubbles  by  breaking  waves,  the  determination  of 
size  distribution,  and  the  temporal  evolution  of  these  features.  Although 
estimates  of  the  volume  fraction  as  a  function  of  space  and  time  are 
possible,  the  direct  measurement  and  knowledge  of  the  volume  fraction  is 
still  lacking.  Laboratory  experiments  with  salt  water  open-ocean  wave 
generation  experiments  as  proposed  by  Longuet-Higgins  would  be  useful  in 
determining  the  single  wave  breaking  plume/cloud  void  fraction  and  acoustic 
characteristics . 
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4.0  NOISE  MECHANISMS 


Breaking  waves  have  been  observed  to  produce  "impact"  noise  as  well  as 
an  ensemble  of  spe. tral  events.  The  spectrum  measured  in  the  deep 
ocean  has  been  shown  to  be  fairly  white  below  -500  Hz  and  in  some 
respects  is  similar  to  the  nondimens ional  spectrum  proposed  by  Franz 
(1959)  to  describe  higher  frequency  sound  generation  by  single  droplet 
impact,  entrained  bubble  oscillations  and  large  scale  wave  breaking 
events.  Prosperetti  (1988)  has  questioned  the  applicability  of  this 
spectrum  at  the  lower  frequencies.  There  is  little  doubt  that  at  the 
higher  frequencies  bubble,  spray  and  splash  are  primary  sources  of 
sound.  However,  at  frequencies  <500  Hz,  individual  bubble  and  splash 
sounds  are  questionable  simply  on  the  basis  of  observed  droplet  and 
bubble  size  distributions. 

Nevertheless,  just  as  a  single  droplet  produces  an  impact  and  a 
subsequent  entrained  bubble  oscillation  at  higher  frequencies,  a 
breaking  wave  may  produce  a  "water  hammer"  type  impact  (Nystuen  (1986), 
Paynter  (1961)),  and  a  subsequent  oscillation  of  the  micro  bubble  plume 
and  cloud  at  lower  frequencies.  This  water  hammer  effect  can  be  shown 
to  produce  a  transient  pressure  pulse  with  a  Fourier  spectral 
amplitude  of  a  dipole,  | P^fR) | -pU^C  CosO/R^  where  p  is  the  density, 

Ujj  the  velocity  spectrum  of  the  impact,  C  the  sonic  velocity,  0 
the  angle  from  the  downward  directed  normal  and  R  the  radial  distance. 
While  wave  breaking  may  not  produce  a  truly  impulsive  event,  the 
resultant  average  spectrum  could  resemble  Franz’s  nondimensional 
spectra.  As  mentioned  by  Nystuen  (1988),  the  velocity  dependence  is 
linear  in  this  case  of  a  water  hammer  impact  as  opposed  to  the  cubic 
dependence  proposed  by  Fitzpatrick  (1959)  for  solid  sphere  impacts. 

4.1  Mixture  Theory 

In  a  bubbly  mixture,  when  the  distribution  of  bubble  sizes  results  in 
resonant  bubble  frequencies  much  greater  than  the  frequencies  of  inter¬ 
est,  then  the  propagation  of  sound  in  the  mixture  can  be  described  by 
an  effective  wave  equation  with  mixture  density  and  sound  speed.  These 
effective  wave  equations  have  been  studied  by  several  investigators 
(i.e.,  Foldy  (1945),  Crighton  (1969),  Van  Wijngaarden  (1968))  and 
confirmed  by  Karplus  (1958)  at  frequencies  <1  kHz.  (See  Carey  (1988).) 

A.  B.  Wood  (1932)  showed  that  at  low  frequencies  in  an  air-water 
mixture  that  when  the  size  of  the  bubbles  and  spacing  between  bubbles 
were  small  compared  to  a  wave  length,  that  for  a  given  volume  fraction 
(x)  the  mixture  density  (pm)  and  compressibility  (*m)  are 
under  equilibrium  conditions  given  by 

Pm  -  (1  -  x)pg.  +  X  P&  and  =  (1  -  x)*a,  +  X 

where  8.  represents  the  liquid  and  g  the  gaseous  component.  The  low 
frequency  limiting  sonic  speed,  C^f ,  is  seen  to  be 

2  dP  ,  -1 

C  -  =  (p  *c  ) 

ml f  dp  mm 
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In  the  case  of  the  air  bubble  water  mixture,  the  process  can  be  con¬ 
sidered  either  adiabatic  or  isothermal.  Since  the  controlling  factor  is 
the  rate  of  heat  transfer  during  bubble  compression  to  the  surrounding 
fluid  which  is  rapid  in  water  due  to  its  large  thermal  capacity,  the 
bubble  oscillations  at  low  frequencies  may  be  considered  isothermal. 

The  following  equations  result: 


„  2^2  22 
-2  _  ji__-  +  um-xA  s  *  P<L  1 

mlf  “  „2  *  2  +  (X)(1  x)  p  C  PflC0 

C0  C  pg  gp8.  l 


-2  -2  -2 
X  -»  0  C  -*  C.  and  x  1  C  -»  C 
ml  mg 


and  when  0 . 002<k<0 . 94 ,  C  ,,  =  P/p.x(l-x). 

mlt  1 


These  equations  show  the  large  effect  on  sonic  speed  as  a  function 
of  void  fraction  and  relative  bulk  modulus  (Fig.  9).  In  the  low 
frequency  range  (<1  kHz)  Karplus  used  an  acoustic  tube  to  determine  the 
standing  wave  pattern  as  a  function  of  air  volume  fraction  and  verified 
these  results.  This  behavior  is  due  to  the  mixture’s  mass  being 
primarily  due  to  the  liquid  while  its  compressibility  is  due  to  the 
gas.  Since  we  are  dealing  with  the  low  frequency  response,  much  less 
than  any  individual  bubble  resonant  frequency,  only  the  volume  fraction 
is  important  and  not  the  bubble  size  distribution. 


Figure  9.  The  relative  bulk  modules  versus  volume  fraction  is  shown 
for  an  air  water  mixture  at  a  pressure  of  one  atmosphere. 
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These  equations  do  not  show  the  dispersive  nature  of  this  sonic 
speed  variation.  If  one  considers  a  bubbly  liquid  with  an  uniform 
bubble  size  whose  corresponding  resonant  frequency  is  o0  and 
dampening  6,  then  one  can  solve  the  component  conservation  equations 
to  obtain  the  following  frequency  dependent  sonic  speed: 

_i  _  u-x)2  _ i _ 

2  2  2  2  2 
C  C  C  (l-u>  /w  +-  2i6u)/u>  ) 

m  P  mlf  o  o 

This  expression  reduces  to  our  previous  equation  for  when  w<<u>0 . 

The  equation  does  show  a  dispersive  character  as  illustrated  in  Fig.  10 
compared  to  the  data  of  Fox  (1955).  Of  interest  is  the  sonic  velocimeter 
range  f>240  kHz,  C->C^ ,  the  sonic  speed  of  the  liquid.  The  dispersive 
character  of  the  phase  velocity  is  the  reason  why  volume  fraction 
measurements  are  necessary  to  determine  the  low  frequency  sound  speed 
changes  near  the  surface  of  the  sea. 
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Figure  10.  The  dispersive  behavior  of  the  sonic  speed  in  air  water 
mixtures  showing  behavior  below,  at,  and  above  bubble  resonance  region. 


4.2  The  Oscillation  of  a  Bubbly  Spherical  Volume 


The  generation  and  scattering  of  sound  from  a  compliant  sphere  immersed 
in  a  fluid  can  be  found  in  classical  texts  on  the  theory  of  acoustics 
(see  Reschevkin  (1963),  or  Morse  (1948)).  In  this  particular  case,  the 
bubbly  sphere  has  no  well-defined  boundary  but  nonetheless  is  localized 
by  perhaps  a  vortex  or  some  other  circulatory  feature  beneath  the 
breaking  wave.  Since  this  analysis  is  somewhat  standard,  only  a  brief 
outline  will  be  presented. 

Here  we  assume  that  the  bubbly  region  (Fig.  11)  is  compact  with  an 
arbitrary  radius  rQ  and  the  region  is  composed  of  micro  bubbles  with 
resonant  frequencies  far  above  the  frequency  of  excitation.  Buoyancy 


18 


forces  and  restoring  forces  such  as  surface  tension  are  not  required. 
The  properties  of  the  bubbly  region  are  described  by  the  mixture  speed 
c  and  density  p  with  a  resulting  compressibility  l/pc^.  The 
micro  bubbles  supply  the  compressibility  and  the  liquid  supplies  the 
inertia.  p 


X->  >r0 


p.  c 


Figure  11.  The  random  collection  of  micro  bubbles  within  a  radii  r 
from  the  origin.  (Pi  =  incoming  wave  or  excitation,  p,c  =  the  pro¬ 
perties  of  the  mixtures,  p, c  =  the  properties  of  water,  and  P  =  the 
radiated  sound.) 

The  source  of  excitation  is  assumed  to  be  global  compared  to  the 
dimensions  of  the  compact  sphere.  This  assumption  means  we  can  assume 
a  plane  wave  expanded  in  terms  of  spherical  harmonics.  The  physical 
reasoning  is  that  the  properties  of  the  bubbly  region  determine  its 
ability  to  radiate  provided  there  is  a  source  of  excitation,  i.e., 

CD 

P.=  P  exp ( iut  -  ikr)  =  P  exp(iut)  l  im(2m+l)  P  (0)J  (kr).,,N 
l  o  o  -  mm  (1) 

m=o 

We  require  continuity  of  velocity  and  pressure  at  the  generalized 
radius,  rQ.  Furthermore  a  radiation  condition  is  imposed  at  large  r  and 
the  field  is  required  to  remain  finite  within  the  bubbly  region. 


The  particle  velocity 

=  (-1/iup)  aPj/ar. 

(2) 

Continuity  of  Pressure 

Pi  +  ps  =  P. 

(3) 

Continuity  of  Velocity 

Vi  +  Vs  =  vn. 

(4) 

The  procedure  is  to  assume  the  scattered  wave  or  radiated  wave  is 
a  sum  of  outward  propagating  spherical  waves. 


P  =  exp (iut)  Ik  P  (0)G  (kr)exp(-ie  kr) .  (5) 

3  m  m  m  m 

m=0 

The  pressure  field  inside  the  volume  is  expanded  in  terms  of 
spherical  Bessel  function  of  the  1st  kind,  i.e., 


P  =  exp ( iut)  l  A  P  (0)J  (kr) .  (6) 

„  mm  m 
m=0 
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When  the  region  is  compact  \  >  2irr0,  we  can  solve  for  the  coefficients 
and  by  using  the  boundary  conditions  (2,  3,  and  4)  and  equating 
each  m  order  terms.  To  O(k-^r^)  ,  we  find  for  Aq  ana  A^  the  following: 


iP  (k  r  )3  (1  -  y/y)/3 

_ Q  Q  Q _  _ 

2  2 
(kr  )  _  (kr  ) 

((1  -  y/y - r — ]  -  i  (kr  )[y/y  - r — ] 

3  O  3 


(7) 


and 


P  (kr  )  (p 
o  o 


p) 


(2 p  +  p) 


(8) 


where  y=pC2  and  y =pC2  .  The  term  Aq  is  interesting  when  the  real  part 
of  the  denominator  is  equal  to  zero,  we  may  say  a  resonance  has  occurred. 


This  occurs  when 


since 


(kr  )2 
o 


2irf 
( - 


r 

o  o 
C 


) 


2 


and 


f 

o 


c2 


yP/p^xd-x)  , 


(9) 


(10) 


which  we  recognize  as  a  modified  Minnaert  formula  for  the  resonant 
frequency  of  a  volume  oscillation  of  a  gas  bubble.  (Note  that  we  have 
inserted  a  factor  y,  the  rate  of  specific  heats.  This  factor  applies 
to  individual  bubble  pulsations  and  for  the  cloud  we  set  y  =  1 
corresponding  to  isothermal  conditions.)  This  result  follows  because  the 
resonant  angular  frequency  is  proportional  to  the  compressibility  of  the 
bubbly  region,  characterized  by  the  stiffness  (4irr0pC2)  and  the  inertia 
(4Trr^p/3).  For  the  case  of  air  with  Ca=,340  m/s  and  pg/pj,  =  0.0013,  we 
find  the  resonant  frequency  of  the  gas  bubble  to  be  f0b  =  6.758  Hz-m/2a 
and  for  the  cloud  with  a  v.f.  of  x  -  0.002  and  Cm  -  200  m/sec,  a 
resonant  frequency  foc=110  Hz-m/2r0.  If  the  average  size  of  each  gas 
bubble  is  50  ym,  then  we  have  f0^=67.6  kHz,  and  if  the  radius  of  the 
air  cloud  is  ro=0.25  m,  we  find  that  foc~220  Hz.  Thus,  the  cloud 
radiates  sound  at  a  resonant  frequency  much  lower  than  its  constituents. 
The  dampening  is  important  to  mention  as  we  expect  dampening  constants 
to  exceed  the  thermal  constant  at  these  frequencies  due  to  the  dirty 
nature  of  sea  water  bubbles  (Filer  (1970)). 

Shown  on  Fig.  12  are  curves  of  resonant  frequency  versus  radius  of 
both  spherical  and  cylindrical  volumes  with  bubble  v.f.  between  10-2 
and  10' ^ .  These  v.f.  are  consistent  with  measurements  of  bubble 
plumes  but  must  be  differentiated  from  the  bubble  densities  and  v.f. 
within  the  near  surface  residual  layer. 
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Figure  12.  The  spherical  and  cylindrical  bubble  cloud  resonant 
frequency  vs  radius  and  void  fraction  (v.f.)  are  shown  and  compared  to 
a  single  spherical  air  bubble. 

The  curves  for  the  spherical  bubble  volume  <»re  taken  from  Eqs.  (9) 
and  (10)  previously  derived  in  this  paper.  Also  shown  are  curves  for  a 
long  cylindrical  volume.  The  estimates  for  the  cylindrical  volume  are 
scaled  by  a  comparison  of  the  resonant  frequency  radius  at  specific 
void  fractions.  This  ratio  of  (Ros/Roc)  the  spherical  radius  to 
cylinder  radius  is  shown  to  be  approximately  3.  The  results  show  that 
the  frequency  radius  dependence  of  a  cylinder  with  a  x  =  0.001 
corresponds  to  a  spherical  cloud  of  x  =  0.01.  These  curves  show  that 
structures  between  10  and  50  cm  can  have  resonant  frequencies  in  our 
range  of  interest  ~100  Hz.  These  calculations  show  that  portions  of 
the  deep  diving  bubble  plumes  have  a  resonant  behavior  within  the 
frequency  100  to  200  Hz  with  reasonable  radial  sizes. 

The  result  of  this  analysis  is  that  micro  bubble  plumes  and  clouds 
to  first  order  can  be  described  by  compact  monopole  sources  of  sound 
since  the  radiation  from  these  bubbly  regions  with  an  arbitrary 
boundary  has  a  natural  frequency  described  by  the  Minnaert  formula. 

4.3  The  Bubble  Plume/Cloud  as  a  Compact  Source 

Minnaert  (1933)  reasoned  that  the  potential  energy  accumulated  during 
the  compression  of  a  gas  bubble  to  its  minimum  radius  is  equal  to  the 
kinetic  energy  of  water  particles  at  the  time  when  half  the  expansion 
has  been  reached,  returning  the  bubble  to  its  equilibrium  radius.  One 
can  perform  this  type  of  calculation  for  either  an  isothermal  or  an 
adiabatic  process.  For  the  individual  gas  bubble,  the  adiabatic  assump¬ 
tion  seems  appropriate.  For  the  oscillation  of  a  cloud  of  bubbles,  one 
should  use  the  isothermal  assumption  due  to  the  thermal  capacity  of 
water . 
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We  have  shown  that  a  bubble  cloud  has  a  first  order  radiation 
characteristic  of  a  pulsating  monopole  source.  Following  the  lead  of 
Strasberg  (1956),  we  represent  this  by  a  simple  harmonic  oscillator 
equation ,  i  .  e  .  , 

2 

V  +  £v  +  WqV  =  0. 

We  assume  that  the  initial  condition  is  V(0)  -=  -vQ  and  that  the  bubble 
cloud  arrived  at  this  condition  in  an  isothermal  process.  The  solution 
to  this  problem  can  be  shown  as 

v(t)  =  -v0  e^d'- 

and  the  radiated  pressure 

po V(t-R/C) 

p(R,t)  =  - -  u(t-R/C). 

From  this  expression,  we  calculate  the  total  energy  in  the  transient 
event  as 


3  2 

-  u  v 

4irC^n  od  o 


where  n  =  £/(»> 


od  ‘ 


If  we  have  Ng  events  per  unit  of  time  and  we  choose  a  reference 
power  corresponding  to  lyPa  at  a  distance  of  1  m,  we  find 


- 


Vi°i 


4ir(lyPa)' 


Using  our  expression  for  the  resonant  frequency  of  the  cloud,  we  find 


E  = 


V/'2" 


1  P  3/2  2 

(  ■)  ("I  v 

3  x  o 


«>/(PR  = 


ST  N_pn  1/2 
4irh 


(4/3irR  ) 
o 


rP ,  3/2  2 

—  [~]  v 

3.X  o 


We  define  a  source  level  (SL)  as  with  representative  values  of  all 
constants  at  lO^pa,  i.e., 


SL  =  10LOG(j>/|pR)  =  10LOG(Nb)  +  20LOG(U  -  15L0G(x)  +  80  dB. 
With  v0  =  ^v0,  ^  =  10^2  and  *~10-4  we  found 
SL  =  10LOG(NB)  +  76  dB . 

These  SL  values  are  large  and  should  be  considered  as  an  upper  bound. 
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This  value  of  the  monopole  source  levels  needs  to  be  examined 
with  respect  to  the  surface  image  interference  effect.  For  a  monopole 
beneath  the  pressure  release  surface,  we  may  show 

| pd ( Ro ) 1 2  =  i P0<R0) I  2  (1  +  V2  +  2v  -  sin2 (kZ  sin©)} 

as  per  the  geometry  shown  in  Fig.  13.  When  the  reflection  coefficient 
is  taken  as  jj  =  -1  and  (Z/X)  =  kZ  <  1,  we  find 


2  2  2  2  2 
|Pd(R0) I  =  |Pq(R0) I  •  16*  (Z/X)  sin©. 


i 


Figure  13.  The  geometry  for  the  surface  image  interference  is  shown 
with  0  taken  as  */\  in  the  downward  direction. 


This  source  expression  is  recognized  as  a  dipole.  Integration  of 
this  over  the  half-space  one  finds 

J>d  =  0>  •  (8ff2/3)  (Z/X)2. 

Our  reference  source  level  now  becomes 

SL  =  10L0G(NB)  +  20LOG( 5 )  -  15L0G(x)  +  20LOG(Z/X)  +  70  dB, 

and  with  Z/\~l/15  we  find  SL~10LOG(NB)  +  67  dB. 

Thus,  for  the  case  of  a  bubbly  sphere  near  the  ocean  surface  with 
5  =  10_2,  a  volume  fraction  of  x~10-4»  and  for  a  frequency  near  100  Hz 
located  at  a  representative  depth  below  the  surface,  sufficient  source 
level  is  achieved.  These  estimates  are  uncertain  as  the  knowledge  of 
the  dynamic  void  fraction  as  a  function  of  depth  and  time  beneath  a 
breaking  wave  is  lacking.  However  our  analysis  shows  that  a  plume/cloud 
can  radiate  sound  with  a  dipole  characteristic  and  is,  thus,  a 
candidate  mechanism  for  the  production  of  sound  at  low  frequencies. 


5.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

This  paper  has  reviewed  data  attained  at  sea  that  demonstrates  that 
breaking  waves  are  a  source  of  low-frequency  (20  to  500  Hz)  sound.  The 
experiments  reviewed  on  individual  breaking  events  showed  short  time 
(5  s)  acoustic  signals  associated  with  breaking  waves.  The  spectra 
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appear  to  be  a  random  collection  of  acoustic  events  in  frequency  and 
time.  The  experiments  reviewed  concerning  deep  ocean  ambient  noise 
measurements  show  a  fairly  white  spectrum  with  a  wind  dependency 
n(200)~1.5.  However,  as  we  have  mentioned  experiments  that  measure 
not  only  noise  but  relevant  oceanographic  and  meteorological  factors 
are  lacking.  We  would  recommend  that  future  ambient  noise  experiments 
measure  not  only  the  10-meter-wind  speed  but  the  air-water  temperature 
difference  and  the  sea  surface  spectrum.  In  addition,  the  sonic 
velocity  variation  near  the  sea  surface  should  be  measured  especially 
as  one  tries  to  characterize  higher  frequency  noise. 

We  have  suggested  that  these  properties  of  ambient  noise  could  be 
explained  by  wave  breaking  impact  noise  (the  water  hammer  effect)  and 
the  collective  oscillation  of  bubble  plumes  and  clouds.  A  key  factor  in 
determining  the  validity  of  this  hypothesis  is  the  time-dependent 
volume  fraction.  Since  no  direct  measurements  are  available  concerning 
this  parameter,  we  must  only  consider  collective  oscillations  as 
possible  noise  sources.  This  hypothesis  could  be  tested  by  examining 
the  sound  scattering  from  micro  bubble  plumes.  We  have  shown  that 
scattering  of  a  plane  wave,  or  excitation  of  a  cloud  of  bubbles  by  a 
coherent  source,  results  in  a  volume  oscillation  that  is  described  by 
Minnaert's  formula  for  the  bubbly  mixture.  Finally,  we  have  attempted 
to  estimate  the  sound  radiation  from  a  cloud  of  bubbles  initially 
compressed  by  the  wave  breaking  impact.  Levels  seem  to  be  adequate  to 
qualify  this  as  a  candidate  mechanism,  but  many  uncertainties  exist  and 
will  most  likely  be  determined  in  future  experiments. 
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I.OW  FREQUENCY  OCEAN  AMBIENT  NOISE:  MEASUREMENTS  AND  THEORY 


William  M.  Carey  and  David  Browning 
Naval  Underwater  Systems  Center 
New  London,  CT,  USA  06320 

ABSTRACT.  Low  frequency  ocean  ambient  noise  data  are  reviewed  and  sum¬ 
marized.  The  experimental  data,  both  omnidirectional  and  directional, 
when  not  dominated  by  shipping  noise,  are  shown  to  suggest  wind  depen¬ 
dent  noise  at  the  low  frequencies  (<500  Hz).  Candidate  mechanisms  are 
examined  with  the  result  that  wave-turbulence  interaction  at  low  sea 
states  and  collective  bubble  oscillations  at  high  sea  states  are  identi¬ 
fied  as  possible  sources  of  this  sound.  A  description  of  the  sonic 
properties  of  bubbly  water  is  presented  for  low  void  fractions  consis¬ 
tent  with  those  observed  in  bubble  clouds  and  plumes  produced  by  break¬ 
ing  waves.  A  description  of  the  collective  bubble-water  mixture  as  the 
resonant  oscillation  of  a  flexible  volume  with  a  sonic  speed  determined 
by  the  properties  of  the  mixture  is  presented. 

INTRODUCTION 

The  interaction  of  the  wind  with  the  ocean  surface  has  long  been  recog¬ 
nized  as  a  major  source  of  acoustic  noise  (Knudsen  (1948),  Wenz  (1962)). 
Measurements  of  the  omnidirectional  noise  at  the  higher  frequencies 
(>200  Hz)  have  been  found  to  exhibit  wind-dependent  characteristics; 
and,  when  not  dominated  by  shipping  noise,  the  most  likely  mechanisms 
are  related  to  bubbles,  spray,  and  splashes  associated  with  white  caps, 
as  well  as  capillary  wave/wave  interactions  (Urick  (1984)).  Furduev 
(1966)  has  proposed  that  the  characteristic  broad  maxima  in  the  ocean 
ambient  noise  spectrum  between  0.2  kHz  and  1  kHz  be  attributed  to  cavi- 
tating  bubbles.  Kerman  (1984)  discusses  these  mechanisms  in  detail 
(also  see  Fitzpatrick  (1959)),  but  stresses  the  noise  generated  by  the 
non-resonant  oscillation  of  entrained  gas  bubbles  which  result  from  wave 
breaking  and  which  are  forced  by  intense  velocity  of  the  gravity-capil¬ 
lary  waves.  For  wind  speeds  with  a  friction  velocity  greater  than  this 
critical  velocity,  Kerman  concludes  that  sound  is  produced  with  a  veloc¬ 
ity  to  the  3/2  power,  frequency  to  the  -2  power,  and  intensity  propor¬ 
tional  to  the  number  of  bubbles.  However,  in  the  absence  of  white  caps, 
since  noise  persists,  capillary  wave/wave  or  non-linear  wave  inter¬ 
actions  may  be  important  (Mellen  (1987),  Kuo  (1968)). 

At  the  other  extreme  of  the  spectrum  (<2-5  Hz),  ambient  noise  asso¬ 
ciated  with  ocean  microseisms  dominates.  Recently,  this  noise  has  been 


1 


LOW  FREQUENCY  OCEAN  AMBIENT  NOISE:  MEASUREMENTS  AND  THEORY 


shown  by  Nichols  (1981)  and  by  Kibblewhite  and  Ewans  (1984)  to  be  due  to 
wave/wave  interaction.  The  microseismic  effect  was  postulated  by  Lon- 
guet-Higgins  (1950)  and  confirmed  by  several  authors,  including  Latham 
and  Nowroozi  (1968).  Several  authors  have  studied  the  generation  of 
sound  through  the  second-order  pressure  effect  ( Brekhovsk ikh  (1967), 
Goncharov  (1970),  Hughes  (1976),  Lloyd  (1981)).  Kibblewhite  and  Evans 
concludes  with  theoretical  arguments  and  measurements  that  the  dominant 
noise  source  in  the  0.1  to  5  Hz  range  is  the  non-linear  wave  interac¬ 
tion.  Although  difficulties  were  found  in  predicting  absolute  levels, 
both  data  and  theory  showed  a  frequency  dependence  to  the  -6  power. 

In  the  very  low  frequency  (VLF,  2-20  Hz)  and  low  frequency  ( LF ,  20- 
200  Hz),  signals  from  surface  shipping  are  a  significant  contributor  to 
the  measured  noise  and  have  been  observed  to  extend  to  500  Hz.  In  this 
region,  noise  contributors  can  be  a  great  distance  from  the  observation 
point,  and  consequently  the  noise  field  exhibits  the  effects  of  sound 
propagation  in  both  the  horizontal  and  vertical  directions  (Carey 
(1986),  Von  Winkle  (1985)).  Wagstaff  (1981)  showed  that,  if  one  knows 
the  locations  and  types  of  ships,  then  one  can  describe  the  characteris¬ 
tic  of  the  horizontal  noise  field.  Although  the  vertical  noise  distri¬ 
bution,  including  the  broad  horizontal  maxima,  could  be  qualitatively 
expla ined,  severa 1  discrepancies  were  observed.  Wind-driven  noise  could 
explain  these  differences,  and  the  sources  of  this  noise  are  the  subject 
of  this  paper. 


EXPERIMENTAL  EVIDENCE 

Omnidirectional  noise  data  at  low  frequencies  which  are  free  from  flow 
and  flow-induced  vibrations  (Strasberg  (1984))  are  very  difficult  to  ob¬ 
tain.  Several  investigators  (figure  1)  measured  the  spectrum  between  2 
Hz  and  2000  Hz  in  the  deep  sound  channel  or  near  the  bottom.  However, 
most  of  this  data  from  the  relatively  heavily  trafficked  northern  hemi¬ 
sphere  reflect  distant  shipping  noise  in  the  2  to  200  Hz  range  and,  con¬ 
sequently,  little  local  wind  speed  dependence  is  observed  such  as  showr 
in  figure  1.  VLF/I.F  ambient  noise  experiments  must  be  carefully  exam¬ 
ined  to  ensure  that  the  results  are  either  from  distant  or  local  sources. 

Wittenborn  (1976)  (figure  2)  performed  an  experiment  with  hydro¬ 
phones  that  spanned  the  water  column.  Hydrophones  within  the  sound 
channel  showed  little  dependence  on  local  wind  speed  between  10  Hz  and 
200  Hz.  However,  the  hydrophone  below  critical  depth  showed  an  inferred 
local  wind  speed  dependence  (10  to  500  Hz)  for  wind  speeds  between  5  and 
15  kns  with  levels  of  47  dB  and  56  dB  re  IpPa  (310  Hz.  The  15  kn  spectra 
showed  a  slowly  vary i n j^t^road  band  characteristic  between  56  dB  at  10  Hz 
and  65  dB  at  500  Hz  (f  ).  Wittenborn  cites  an  earlier  experiment  with 
noise  levels  of  69  dB  for  300  Hz  at  30  kns,  compared  to  the  300  Hz  lev¬ 
els  of  44  dB  at  5  kns,  51  dB  at  10  kns,  and  63  dB  at  15  kns.  These  re¬ 
sults  suggest  two  wind  noise  mechanisms  for  the  cases  of  low  and  high^ 
sea  states  with  the  intensity  having  a  squa red^vgloc ity  dependence  (U  ). 
The  abrupt  transition  between  10  and  15  kns  (U  '  ,  based  on  the  levels 
at  10  and  15  kns  and  not  considered  a  true  velocity  dependence),  as 
shown  in  figure  2,  may  be  a  threshold  characteristic  associated  with 
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Figure  1.  Selected  Low  Frequency  Ambient  Noise  Measurements.  The  re¬ 
gion  less  than  5  Hz  is  dominated  by  wave/wave  interaction.  The  measure¬ 
ments  between  5  Hz  and  300  Hz  show  little  local  wind  speed  dependence, 
but,  rather,  the  effects  of  distant  shipping  and  other  distant  noise 
sources . 


Figure  2.  Ambient  Noise  Level  vs.  Frequency  for  the  Wittenborn  Experi¬ 
ment.  The  4850  m  deep  hydrophone  shows  the  local  wind  speed  dependence 
(=200  Hz)  with  the  influence  of  distant  noise  sources  less  than  100  Hz. 
The  3960  m  deep  hydrophone  is  dominated  by  distant  noise  sources  less 
100  Hz. 


breaking  waves.  (These  results  agree  with  the  observations  of  Worley 
(1982),  insofar  as  his  data  show  a  threshold-type  behavior  between  the 
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NL  =  NL  i  +  20nLOG(W  S  ) 


INVESTIGATOR 

— 

f 

n 

w  s  (m 'sec.) 

PIGGOTT  JASA  36(11) 

13 

2  1 

10-20 

13 

0 

1-5 

141 

1  53 

3  5-20 

<  50 

2  1 

PAYNE  JASA  (1967) 

50-100 

2  2 

WHITTENBORN  0976) 

200-300 

1.65 

2  5-5 

34 

5-7  5 

1 

7  5-15 

CROUCH  JASA  3(2)-72 

11 

1.73 

5-20 

28 

2  09 

SHOOTER  JASA  73(6)-8l 

1-1  39 

5-10 

1  1-1.32 

10-15 

BB 

1  36/1.57/  81 

WORLEY  JASA  71  <4)-82 

100 

85-1  5 

10-15 

200 

1  65-2  0 

5-10 

WILSON  JASA  7301-83 

10 

2  07 

5-10 

BURGESS  JASA  7301-83 

37 

1  66 

5-15 

FACTORS  AFFECTING  WIND  SPEED  DEPENDENCE  ARE  1  DISTANT  SOURCES. 
2  OVERLAPPING  WIND  SPEED  REGIONS.  3  SOUND  PROPAGATION  FACTORS 


Table  I.  Low  Frequency  Ambient  Noise  Wind  Speed  Dependence 


data  corresponding  to  wind  speeds  between  2.5  and  5  kns  and  between  5 
and  10  kns  at  200  Hz.  This  effect  was  especially  pronounced  at  400  Hz.) 

Although  Wittenborn  made  use  of  both  refractive  effects  and  bathy¬ 
metric  blockage,  noise  from  distant  sources  was  still  found  to  influence 
his  results  (for  example,  see  figure  2  between  10  and  100  Hz).  The  cor¬ 
rupting  influence  of  distant  noise  sources  (ships,  whales,  volcanoes, 
etc.)  has  the  effect  of  obscuring  the  low-frequency  local  wind  speed  de¬ 
pendence.  Consequently,  the  literature  reveals  a  variety  of  estimated 
wind  speed  dependencies;  i.e.,  the  estimate  of  a  parameter  n,  where 
NL  =  NLj  +  20  n  log  (W.S.).  (The  mean  square  pressure  would  increase 
with  2n  power  of  wind  speed.)  Table  1  lists  several  of  these  estimates 
of  n,  ranging  from  0.85  to  2.0  for  wind  speeds  between  10  and  20  m/sec. 
The  problem  with  these  estimates  also  lies  in  the  fact  that  the  data 
clearly  show  a  region  of  no  wind  dependence,  a  threshold-type  behavior, 
and  region  with  a  wind  dependence  of  n-2.0. 

Figure  3  illustrates  this  trend  with  the  data  of  Piggott  (1964). 

One  observes  the  frequency  dependent  cross-over  between  the  low  wind 
speed  and  higher  winds  regions.  Furthermore,  the  lower  the  frequency, 
the  higher  the  wind  speed  will  be  at  which  the  wind  speed  dependence 
point  is  observed. 

Distant  noise  sources  influence  vertical  noise  directionality  (Von 
Winkle  (1985),  Browning  (1982),  Bannister  (1986)).  This  influence  of 
the  distant  source  produces  a  broad  maximum  in  the  vertical  noise  inten¬ 
sity  centered  on  the  horizontal.  This  phenomenon  results  from  the  con- 
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Figure  3.  Ambient  Noise  Spectrum  Level  vs.  Wind  Speed,  Piggott  (1964) 


version  of  higher  angle  rays  to  lower  angle  rays  by  either  reflections 
from  the  basin  boundaries  and  seamounts  or  refractive  effects  due  to 
shallowing  sound  channels  at  the  higher  latitudes.  Wagstaff  attributed 
this  effect  to  surface  ships.  Since  the  spectral  variation  of  the  hori¬ 
zontal  noise  is  generally  smooth,  and  since  ship  signatures  are  narrow- 
band  in  this  frequency  range,  wind-produced  noise  over  seamounts, 
slopes,  and  at  high  latitudes  was  speculated  to  be  an  important  contri¬ 
butor.  The  broad  maximum  along  the  horizontal  has  been  observed  in  var¬ 
ied  geographical  locations,  such  as  the  sparsely  shipped  Southern  Hemi¬ 
sphere  waters  of  the  South  Fiji  Basin  (shown  in  figure  4).  At  the  lower 

frequency  the  data  clearly  show  a  broad  maximum.  At  105  Hz  one  observes 

the  influence  of  a  single  ship.  These  results  are  similar  to  data  ob¬ 
tained  in  the  North  Pacific  and  the  North  Atlantic  (Carey  (1986)). 

The  experimental  data  were  examined  to  obtain  measured  levels  use¬ 
ful  in  the  estimation  of  the  source  level  of  wind-produced  noise  at  the 
sea  surface.  These  results  are  shown  in  Table  II,  primarily  at  50  Hz. 

The  estimated  levels  based  in  the  Wittenborn  data  are  shown  in  the  table 

to  be  between  43  dB  at  5  kns  and  51  dB  at  15  kns,  consistent  with  the 
estimates  by  Wilson  and  Kewley  using  the  same  data.  Vertical  noise  can¬ 
not  be  used  for  local  wind-driven  noise;  however,  estimates  for  a  cylin¬ 
drical  basin  with  sloping  sides  yields  levels  in  the  50  dB  range.  Kew¬ 
ley  has  carefully  estimated  source  levels,  and  his  curves  are  shown  in 
figure  5. 

In  summary,  we  have  presented  data  which  indicate  the  presence  of  a 
wind-driven  noise  in  the  10  to  200  Hz  region  of  the  spectrum.  The  low 
wind  speed  range  (<^8-10  ra/sec)  appears  to  have  a  weak  dependence  on  the 
wind  speed ,  0  <  n  <  1;  the  high  wind  speed  region  (>7.5  to  15  m/sec) 
appears  to  have  a  dependence  of  0.85  <  n  <  2.  These  estimates  point  to 
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Figure  4.  Vertical  Noise  Spectrum  Level  versus  Angle  from  the  Horizon¬ 
tal,  Browning  (1986) 


the  uncertainty  in  our  knowledge  of  wind  speed  dependence  and  spectral 
characteristics . 


POSSIBLE  MECHANISMS 

The  fundamental  mechanisms  for  the  production  of  sound  in  turbulent 
regions  may  be  derived  from  first  principles.  The  basic  procedure  can  be 
found  in  several  treatments  on  hydrodynamic  noise,  most  notably  Light- 
hill  (1979),  Ffowcs  Williams  (1969),  Dowling  (1983),  and  Ross  (1976). 

We  have  rederived  the  inhomogeneous  wave  equation  with  source  terms  in 
appendix  A  for  the  purpose  of  ranking  the  various  mechanisms  capable  of 
the  production  of  sound  at  the  surface  of  the  sea  in  the  10  to  200  Hz 
range.  The  basic  approach  is  to  write  ^he  equations  governing  the  con¬ 
servation  of  mass  and  momentum  with  source  terms.  The  equation  of  state 
is  specified,  fluctuation  quantities  assumed,  linearization  employed, 
and  the  inhomogeneous  wave  equation  is  formed.  The  integral  solutions 
to  this  equation  are  formulated  by  use  of  the  Kirchoff  Method  (Stratton 
(1941),  Jackson  (1962))  and  of  the  divergence  theorem.  The  derivation 
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•  OMNIDIRECTIONAL  MEASUREMENTS  WITH  HYDROPHONE  BELOW  CRITICAL  DEPTH 


INVESTIOATOR/FREQ. 

10  Hz 

50  Hz 

100  Hz 

w.s.(kns) 

WHITTENBORNE  (1982) 

50 

44 

5 

N.E  PACIFIC 

5  i 

47 

10 

58 

60 

15 

MORRIS  (1978) 

■bH 

70 

63 

10 

N.E.  PACIFIC 

•  VERTICAL  NOISE  MEASUREMENTS  ALONG  THE  HORIZONTAL: 


INVESTIGATOR/FREQ. 

50  Hz 

70-80  Hz 

100  Hz 

BROWNING  (1982) 

67  dB 

70 

65 

FIJI  BASIN 

WAGSTAFF  (1981) 

— 

— 

69-66 

N.E.  PACIFIC 

WALES  (1981) 

95 

87 

82 

N.W.  ATLANTIC 

AXELROD  (1965) 

— 

65-69 

N.W.  ATLANTIC 

FOX  (1964) 

— 

— 

60-65 

N.W.  ATLANTIC 

FISHER  (1986) 

74 

— 

60 

E.  PACIFIC 

ANDERSON  (1979) 

- 

- 

65-69 

N.  PACIFIC 

•  SEMI-EMPIRICAL  ESTIMATES:  * dB//(uPa2/Hz ) 


WILSON  (1983)  NL=50dB+30LOG(W.S.)  @  50  Hz 

BURGESS  &  KEWLEY  (1983)  Nl.=37.3  +15LOG(W  S.)  @  50  Hz 

•  SOURCE  LEVEL  ESTIMATES  BASED  ON  OMNI  MEASUREMENTS: 

1=  10  Hz  41  (dB//uPa*<2/Hz  @  1m)@  5  kns;  50  dB  (3>  15  kns 
50  Hz  43  dB  (3)  5  kns  51  dB  <®  15  kns 

100  Hz  37  dB  @  5  kns  56dB@15kns 

•  VERTICAL  LEVELS  YIELD  FOR  A  MEAN  WIND  SPEED  OF  lOkns,  50dB@50Hz  AND  54-56dB@100Hz 

Table  II.  Ambient  Noi6e  Source  Levels 


In  appendix  A  If.  similar  to  those  of  Huon-I.I  (1981)  and  Yen  (  1979),  anti 
the  basic  result  is  the  following: 

AttCo(P-Po)  =  4ttP  =  [  |0q/3t]dV/R  -  0/c'x.  |  [F,]  dV/R  +  — ^  /  [T„]  dV/R  - 

-  I  l,[9pU,/ut]  dS  +  0/9x,  |  l,|2p0U,u,  +  p0u,u,  +  P,,]  — 4~—  • 

R 

The  first  term,  /  [rjq/clt}  dV/R  ,  represents  a  monopole  term,  q  represents 
mass  addition  rate  per  unit  volume.  The  second  term  represents  an  ex¬ 
ternal  force  acting  on  the  volume  and  has  a  dipole  character.  These  two 
terms  could  be  important  in  the  incorporation  of  entrained  bubble  oscil¬ 
lation  and  translation.  The  third  term  is  the  Lighthill  turbulence 
stress  tensor  and  Is  known  to  represent  an  acoustic  quadrupolc.  The 
term  f  l,(r)/yU./0t]  dS  Involves  the  motion  of  the  boundary  and  can  act  as  a 
monopole.  The  final  integral  involves  the  turbulent  and  compressive 
stresses  acting  on  tie  boundary  and  is  seen  to  have  a  dipole  character. 
In  particular  the  term  d/3x,  /  l,[2poU,u,]  dS/R  represents  the  wave  turbulence 
Interaction  and  is  dominant  since  it  represents  a  nroduct  of  a  first 
order  U,  ;  and  second  order  term  u,  . 

NoIbc  generation  by  the  interaction  of  surface  waves  and  turbulence 
near  the  surface  wan  suggested  by  Goncharov  (1970).  calculated  lev¬ 
els  of  80  dB  at  10  Hz  and  40  dB  at  100  Hz  by  assuming  a  Pierson-Mosko- 
witz  surface  wave  spectrum  and  Kolmogorov's  similarity  hypothesis.  His 
expression  can  be  shown  to  be  equivalent  to  the  above  integral.  Howev¬ 
er,  instead  of  using  velocities,  be  employs  the  displacement  spectrum 
for  the  surface  wave  and  turbulence.  His  expression  is  p(w)  A0rt?/rjp . 
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Figure  5.  Source  Level  vs.  Frequency  and  Wind  Speed,  Kewley  (1986) 


Yen  and  Perrone  (1979)  derived  expressions  yielding  the  frequency-depen¬ 
dent  radiation  characterlst i cn  for  the  wave/wave,  wind/turbulcncc,  and 
wavc/turbulonce  interaction  mechanisms.  Their  results  for  the  wave/tur¬ 
bulence  interaction  (70  dB  at  10  Hz  and  50  dB  at  100  Hz)  show  a  linear 
dependence  on  surface  wave  velocity  (U)  and  an  inverse  square  dependence 
on  frequency  (to  7)  : 

P(w.k)  -  2  •  10' Vcos40  •  U /<o‘  ■ 

The  Yen  &  Perrone  result  contains  three  interesting  factors.  The 
linear  dependence  on  surface  wave  velocity  is  consistent  with  the  pre¬ 
viously  discussed  experimental  results  prior  to  wave  breaking.  The  w  7 
dependence  is  also  consistent  with  the  observed  behavior  at  low  frequen¬ 
cies;  i.e.,  an  overlap  region  composed  of  the  interaction  of  the  low  and 
higher  frequency  roll-offs.  However,  of  particular  note  is  the  COsV  de¬ 
pendence.  This  sharply  peaked  angular  dependence  would  accentuate  the 
role  of  the  ocean  bottom  and  basin  boundaries  with  respect  to  the  verti¬ 
cal  noise  directionality.  Thus,  wave/ turbulence  interaction  could  be  a 
source  of  noise  in  the  2  to  200  Hz  region  for  those  sea  states  low 
enough  that  breaking  waves  do  not  occur,  due  to  the  fact  that  it  appears 
ns  a  physically  realizable  mechanism  (considering  the  uncertainty  of  the 
turbulence  spectrum). 

Kerman  (1986)  shows  that,  above  a  critical  wind  speed  of  approxi¬ 
mately  10  m/s,  small  (micrometer  (pm))  bubbles  arc  produced  and  can  be  a 
source  of  higher  frequency  sound.  Thorpe  (1986,  1982)  has  performed  in¬ 
teresting  experiments  which  demonstrate  the  existence  of  bubble  plumes 
and  layers  composed  of  pm- size  Jjubbles^fmonn  bubble  size  approximately 
50  pm  with  densities  between  10  to  10  bubbles/m  )  extending  several 
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meters  below  the  surface.  M.Y.  Su  (1984)  has  shown  that  fresh  water 
breaking  waves  produced  in  wave  tanks  produce  bubble  plumes  which  pene¬ 
trate  to  depths  on  the  order  of  signicant  wave  height,  with  bubbles  of 
centimeter  diameters  due  to  ccalesence  (which  is  absent  in  salt  water). 
Several  reviews  (see  bubble  refeiences)  have  been  written  on  the  exis¬ 
tence  and  densities  of  bubbles  produced  by  breaking  waves.  At  high  sea 
states  a  residual  layer  is  formed  of  micro-bubbles  with  a  density  that 
decreases  exponentially  with  depth,  and  bubble  plumes  which  are  con- 
vected  to  several  meters  depth  by  the  vorticity  beneath  the  wave.  Even 
though  individual  bubble  oscillations  with  these  micron-size  bubbles 
could  not  produce  VLF/LF  noise,  collective  oscillations  of  the  bubbly 
mixture  driven  by  the  hydrodynamic  pressure  field  could  produce  sound. 

It  is  well  known  (see  appendix  B)  that  a  small  amount  of  bubbles  in 
water  significantly  changes  the  bulk  compressibility  while  not  drasti¬ 
cally  changing  the  density.  Wood  has  shown  that  the  sonic  velocity  (Cm) 
can  be  described  by  the  following  relationship  between  void  friction 
(x),  density  (pm),  and  bulk  compre s s ib i 1 i te s  (K): 

C,/  =  ((1 -x)  p,  +  xpg]((1-x)  K  +  xK,;]  ■ 

The  consequence  of  this  result  is  shown  in  the  figures  of  appendix  B. 
Small  volume  fractions  result  in  large  changes  in  the  sonic  speed  when 
the  mixture  can  be  treated  as  a  continuum.  For  example,  the  sonic  ve¬ 
locity  of  the  bubble  mixture  with  a  0.27,  volume  fraction  is  approxi¬ 
mately  225  m/s.  Ffowcs  Williams  (1969)  describes  the  efficiency  of  the 
radiation  from  a  cloud  of  bubbly  turbulent  flow: 

4  nCa(p~pQ)  =  4  rr P  =  I  (dq-'3t)  dV.  R  • 

For  a  compact  source  with  a  small  gas  volume  fraction 

AnCo(p-po)  -  477P  i  q  dV, 

R  c\ 

Ffowcs  Williams  estimates  that 

q  =  -p  ! n ( 1  - x )  —  -Cn''  DP  Dt 

dp  =  Cm  A(1-x)p  =  -pCm  Ax  -  p-po  [p/ R)  m4  (Co/Cm)?.  m  =  u  'Co 

Thus  he  concludes  that  "a  cloud  of  bubbly  flow  radiates  very  much  more 
efficiently  than  turbulence  alone;"  that  is,  the  radiation  from  such  a 
flow  would  be  (C/Cm)  -  1975-times  larger  than  the  radiation  from  turbu¬ 
lent  flow.  However,  one  must  account  for  the  presence  of  a  pressure 
release  surface. 

An  alternative  approach  is  to  consider  the  bubble  cloud  as  a  flexi¬ 
ble  sphere  of  radius  a  with  composite  mixture  properties  and  to  assume 
it  is  compact  with  respect  to  the  acoustic  wave  length  and  the  vorticity 
and  turbulence  scales.  Then  the  forced  oscillation  of  the  bubble  cloud 
in  absence  of  a  boundary  is  t' =  t-R/Co  »  Q  =  |  q  dV  =  ^  (p0v) 

dt 

V(t)  3o/a  3p  , 

P(R.t)  =  Q(t')/47rR  =  p0  ---  = - - - nr  (Co/Cm)  f(w)/R, 

AttR  (1-x) 

where  f(w)  represents  the  simple  harmonic  oscillator  transfer  function. 
Thi6  forced  oscillation  of  a  bubble  cloud  can  have  a  resonant  behavior 
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(Fitzgerald  and  Mellen  (1982)).  The  bubble  cloud  is  simply  a  monopole 
source  and  the  presence  of  the  boundary  can  be  approximately  taken  into 
account  via  the  surface  image  interference  effect.  Thus,  we  find  the 


following  for  a  cloud 
face  driven  in  forced 

Ip(R.t)!2  ssf"- 


of  micro-bubbles  below  the  pressure  release  sur- 
oscillation  by  the  hydrodynamic  forces: 

2  4  (  2ttz 
R5 


3cj?a3p 


(1-x) 


m^Co/Cm)"  f(w 


ui  v/ujruauiAv  aw 


sin20 


This  expression  shows  a  frequency-dependent  efficiency  approximately 
(z/X)  at  a  given  bubble  cloud  depth.  This  term  indicates  that,  as  seas 
pick  up,  the  deeper  the  plume,  the  more  efficient  the  radiation  at 
longer  X.  Furthermore,  we  note  that  this  monopole  has  an  m  improvement 
over  the  non-compact  bubble  cloud. 

Thus,  the  low-frequency  noise  could  be  caused  by  wave  turbulence 
prior  to  wave-breaking  and,  thereafter,  by  aggregate  bubble  (bubble 
cloud)  oscillations  exhibiting  a  threshold-type  behavior  and  velocity 
squared  dependence  for  the  mean  radiated  pressure. 


APPENDIX  A:  DERIVATION  OF  THE  SOURCE  INTEGRALS 

The  purpose  of  this  appendix  is  to  briefly  outline  the  derivation  of 
source  terms  important  to  the  production  of  sound  near  the  surface  of 
the  sea. 

Conservation  of  mass :  dp'/d\  +  0p'v,/dx,  =  q  • 

.  c  .  dp'v!  0p’v,v,  dP'„ 

Conservation  of  momentum:. JL —  +  — _ — —  = - !_  +  |-0,  , 

0t  dx,  3x, 

where  P’1(  =  -p'<5,,  +  pD,,  +  p,0<5,,  (ref.  Hinze,  p.  17) 

p,  =  2/3 p,  D,,  =  0U,/0x,  +  dll./dx,  and  0  =  1/2D„  =  9Uu/3xh  • 

Taking  9/3t  of  the  continuity  equation  and  3/9x,  of  the  momentum  equation 
yields  upon  subtraction: 

2  2  ?2p'v,v |  02P',,  9Fei  ^  9q 

^  ^  3x,  3x,  dx,  0 xf  0x,  3t 

dp o  dpo 

v  =  U  +  u,  (0U/0x.  =  0),  p  =  po  +  P.  — — =  -  0 

d?p/d f  -  Co  dip!dx\  -  U.U,  d’p/dx,  dx,  + 

d?(po+p)n.u,  _ 29?(po+p)U,u, 

dx,  dx,  3x,  0x, 

dClp6,,  0P.|  3F„  dq 

-  ....  4-  - -  •  . .  4  ■  r-  — 

0x,  0x,  3x,  0x,  3x,  d\ 
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For  the  case  of  incompressible,  invisid  flow  with  no  sources  or  sinks: 

y  -  „  C2  Po  U,U,  _j  . 

P  pou.u./'x  ex  +  2 - +  £  pi  ex  .  =  0  ; 

ox.  ex. 

Compressible  fluid: 

pCl  =  P 

1/C2  p2P'3t2  -  02P/  Jx2  -  c:q it  -  ■;■<,.  'ex  +  "'T  /Jx  x. 

T,  =  2p0  U  u,  +  pou.u,  +  P„  -  Cip6,  . 


(in  most  instances  P,p  -  CoP<5,,  =  0  ,  (Lighthil  1 ) ) . 

Since  P  =  po  3<///Jt  =  -\u>  pr  <//  ,  we  finally  have  the  wave  equation  with  the 
source  terms: 

-3  1  d2ip  1  \  cq  3fe,  32T.,  / 

- 5 —  -  -  -~r-  — j-  = - (  — —  -  — 7 - +  — — : —  (  =  -477  f(x.t)  . 

ex  Co  't  icvpo  |  et  ex  ex  ex,  \ 


This  inhomogeneous  wave  equation  can  be  integrated  by  use  of  the 
Kirchoff  method  (Stratton  (1941)  and  Jackson  (1962))  to  yield: 

I  dV  [f]  1 

<//( x . t)  -  - +— —  I dS  [1/R  dili/ or  -  3/3n(1/R )ip  +  1/CoR  3R''3n  dtl//d\] 

x  -  x'  4tt 

This  solution,  when  applied  to  cur  specific  problem  with  the  properties 
of  3[  ]/3x^  and  3[  ] / 3 y ^ ,  as  well  as  the  divergence  theorem,  yields  the 
desired  results: 

4rrCo(p  -  Po)  =  4rrP  =  J  [dq/'d\]  dV/R  -  dlex.  |  |U,)  dV/R  +  32/3x,3x,  j  IT,]  dV  'R  - 

-  I  I,  [3pu,/3t]  dS  +  3/3 x.  I  I,  [2p0U,u,  +  pou,u,  +  P,,]  dS/R  . 

APPENDIX  B:  MIXTURE  THEORY 

A. B. Wood  (1932)  showed  that  the  sonic  speed  could  be  calculated  for 
an  air-bubble/water  mixture  by  use  of  the  mixture  density  (pm)  and  the 
mean  compressibility  (km).  The  mixture  can  be  treated  as  a  continuous 
medium  when  the  bubble  diameter  (d)  and  spacing  between  the  bubbles  (D) 
are  much  less  than  the  wavelength  of  sound.  In  the  case  of  low  frequen¬ 
cies,  for  the  mixture  with  a  volume  fraction  (X)  of  gas  we  can  calculate 
the  mean  density  and  compressibility  as  follows: 

Pr n  =  ;i  '  X)p.  +  XPq 

-dVm  dVi  V,  d  Vq  V9 

K  (n  =  — —  = - -  pr  —  +  — - - =  (1  -  x)Ki  +  x  Kg  . 

V,r  dP  V:  dP  Vm  v,dP  vm 

This  imples  that  a  state  of  equilibrium  prevails  and  the  mixture  mass  is 
conserved,  and  the  pressure,  P,  is  uniform  throughout  the  mixture  (a  low 
frequency  assumption).  Since  the  sonic  speed  is 

C2  -1  dP/dp  =  (pK)  ’  , 


we  have 


C  n  =  Cmll  =  ((1  -  X)p,  +  XPg]  [(1  -  X  )  K  +  XKg 


Crr2  -  (1  -  X)/C?  +  X2/Cg  +  ( X)  (  1  -  X) 


Pq  C2  ♦  P?C- 
pi  Pg  C  Cg 
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Figure  B-l.  Low  Frequency  Approximation  for  the  Mixture  Speed 


The  expression  for  the  sonic  speed  poses  the  question  of  whether 
the  gas  compressibility  is  described  by  an  isothermal  or  adiabatic  pro¬ 
cess,  especially  since  the  single  phase  sonic  speed  is  known  to  be  adia¬ 
batic.  However,  in  the  case  of  an  air-bubble/water  mixture,  the  con¬ 
trolling  physical  factor  is  the  transfer  of  the  heat  generated  in  bubble 
compression  to  the  surrounding  liquid.  If  the  transfer  is  rapid, then 
the  bubble  oscillation  is  i sotherma  1 ,  Mv/r)P  =  -v/P  KQ  =  1/P)  ,  as  compared  to 
the  adiabatic  condition  (dv/dP  =  -v/yP.  Kg  =  1/yP)  •  Thus,  in  use  of  the 
above  equations  one  must  use  either  for  the  adiabatic  or  isothermal 
case,  C9,  =  Cga/v  Y  •  Isothermal  conditions  are  most  likely  to  prevail  for 
air-bubble/water  mixtures  due  to  the  large  thermal  capacity  of  water. 
Examination  of  the  above  expressions  shows  that  as  x  -  0,  Cm2  -  Ci  2,  and  as 
x  —  1 ,  Crp2  —  Cq2>  as  one  would  expect.  The  striking  characteristic  revealed 
by  these  equations  (shown  in  figure  B-l)  is  the  sharp  reduction  in  the 
sonic  velocity  at  small  volume  fractions;  i.e.,  X“0.002  -+■  ^"225  m/sec. 
These  equations  may  be  approximated  for  the  air/water  mixture: 

ca  _  Xp  x  _ P 

p.x  n  -  X)  p,x  (1  -  X) 

C.T.  (x  =  0  5)  =  20  m/sec 

Karplus  (1958)  used  an  acoustic  tube  to  determine  the  standing  wave  pat¬ 
tern  as  a  function  of  air  volume  fraction.  His  results  are  shown  in 
figure  B-2.  Close  agreement  was  found  between  the  inferred  sonic  speeds 
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PERCENT  AIR  BY  VOLUME  FREQUENCY 

VELOCITY  OF  SOUND  VELOCITY  OF  SOUND 

IN  WATER  CONTAINING  AIR  BUBBLES  in  WATER  CONTAINING  AIR  BUBBLES 

2500 


2000 

1500 

Cp»(m/s) 

1000 


500 


0  30  60  90  120  150  180  210  240  270  300  330  360  390 

f  (kc/s) 

SINGLE  RUN  PLOT  OF  PHASE  VELOCITY  vs  FREQUENCY 


Figure  B-2.  Measured  and  Computed  Mixture  Sonic  Speeds.  a)  Measured 
Mixture  Speed  vs.  Volume  Fraction  at  500  and  1000  Hz,  Karplus  (1958);  b) 
Measured  Dispersive  Character  of  Mixture  Speed  Below  2000  Hz,  Karplus 
(1958);  c)  Measured  and  Computed  Mixture  Sonic  Speed  Showing  the  Beha¬ 
vior  Below,  At,  and  Above  Resonance,  Fox,  et  al.  (1955) 
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and  results  calculated  with  Wood's  expressions.  Similar  results  have 
also  been  observed  at  the  low  frequencies  by  Campbell  and  Pitcher 
(1955).  These  results  are  also  observed  at  the  higher  frequencies  above 
and  below  resonance.  Several  studies  and  texts  have  been  written  on 
this  subject  and  are  listed  in  the  references.  An  example  of  the  agree¬ 
ment  between  theory  and  measurement  near  the  vicinity  of  bubble  reso¬ 
nance  is  shown  in  figure  B-2c.  It  is  important  to  note  that  most  calcu¬ 
lations  performed  at  these  higher  frequencies  use  Km  =  K^  +  Kg,  rather 
than  the  Wood  approach  Km  =  ( 1  -  x)K^  +  xKg.  This  difference  is  unim¬ 
portant  near  resonance  and  for  small  volume  fraction  but  is  important 
as  one  approaches  the  low  frequencies  of  interest  to  this  paper.  One 
can  show  that  the  correct  expression  is: 

_  n  -  xr  _  i 

C?  C'  Ci'i  ((1  -  Co' '  wcl  +  2io  (jj  ll><.  i 

when  h.f.  and  l.f.  are  the  high  frequency  and  low  frequency  values  of 
the  sonic  speed. 
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